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Abstract. We consider nonlinear parabolic equations of the type 

u t - div a{x,t,Du) = f(x,t) on(] T = Ox (~T,0), 

under standard growth conditions on a, with / only assumed to be inte- 
grable. We prove general decay estimates up to the boundary for level sets 
of the solutions u and the gradient Du which imply very general estimates 
in Lebesgue and Lorentz spaces. Assuming only that the involved do- 
mains satisfy a mild exterior capacity density condition, we provide global 
regularity results. 
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1. Introduction 

The aim of the present paper is to establish global regularity results for the 
solutions of the following class of Cauchy-Dirichlet problem 



XI) 



ut — dwa(x, t, Du) = f(x, t) in Qt = ^ x (— T, 0), 

U = on dpar^T, 



where £1 is a bounded open set in H n , n > 2, T > 0, <9 pa r^T is the usual parabolic 
boundary of £It- Here we assume that /, in the most general case, is a summable 
function in Sly, and a : £It x R n — > R n is a Caratheodory regular vector field 
satisfying standard monotonicity and growth conditions; i.e., 



1.2 



v\z2 - zi | 2 < (a(x, t, z 2 ) - a(x, t,z\),z 2 - z\) 
\a(x, t,z)\ < L\z\ 



for every z±, z 2 , z G M n and (x,t) G Ut] the structure constants satisfy < v < 
1 < L. 

We will focus here on the case when / belongs to the Lebesgue space L 7 (r2^) 
in a range of 7 that does not necessarily permit to obtain the existence of finite 
energy solutions u 6 L 2 (— T, 0; Wq' 2 (Q,)) to problem (1.1). 

Remark 1.1. In this paper we are not going to consider any solution to the problem 
(1.1), but rather we shall deal with the (very) weak solutions u € L l {—T, 0; Wq 1 ' 1 ^)) 
obtained via the Boccardo-Gallouet by now standard approximation procedure ([9]; 
see forthcoming Section 2.1). These solutions are called SOLA (Solutions Obtained 
by Limits of Approximations) and turn out to be unique in their class in the case 
/ € L 1 (f2'r), see [12, Theorem 4.1]. For such reasons, with some abuse of termi- 
nology, we shall call in the following the unique SOLA to (1.1) "the solution" to 
problem (1.1). 

In the classical Lebesgue framework, Boccardo, DalPAglio, Gallouet and Orsina 
(see [8, Theorem 1.9]) show the existence of a unique solution 

(1.3) u € L 9 (-T, 0; W£' q (n)) with q = j^- 

to problem (1.1) under the assumption (1.2), provided that the datum / satisfies 

/ G L 7 (r2y) for some 1 < 7 < 



N + 2 

Here ./V := n + 2 denotes the homogeneous parabolic dimension. Moreover, the 
solution u belongs to L a (Qx) with a = Nj/(N — 2-y). This result is optimal in 
the scale of Lebesgue spaces; however it is natural to ask for a more accurate 
scale to describe the regularity of Du in dependence on the datum /. In this 
respect, recently, the results above were extended in all the most familiar spaces 
of rearrangement (Lebesgue, Lorentz, Orlicz) and non-rearrangement once, up to 
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Lorentz-Morrey 1 , by Habermann and the first author of this paper in [5] (see, also, 
[25, 14, 26]). We would like to remark that a fine analysis in Lorentz spaces regular- 
ity is often decisive to establish higher regularity, as it has been shown, for instance, 
by recent developments in the theory of harmonic and biharmonic mappings (see, 
e.g., [41, 22, 35] and the references therein). Precisely, in [5, Theorem 6.6], starting 
from the datum / in the Lorentz space L(7,<7)(J1t), for 1 < 7 < 2N/(N + 2), the 
authors prove that for any < q < 00 



extending to the parabolic setting a very general approach introduced by Min- 
gione in [31, 32, 33] for the regularity of solutions to the elliptic counterpart of 
problem (1.1). 

After [31], many important extensions and generalizations have been considered 
in the literature, in which a similar approach based on maximal operators tech- 
niques has been pursued to achieve regularity results in nonlinear degenerate prob- 
lems: see, e.g., [33, 13, 15] for elliptic problems, [1, 5, 14, 17, 26, 27] for parabolic 
ones and the references therein. It is important to stress that the aforementioned 
papers, despite considering more and more general settings, provide local results 
not considering the regularity properties up to the boundary of the solutions. 

Recently, some attention has been devoted to such global extension, as in the 
paper by Phuc [40], where Lorentz regularity in the whole domain Q has been 
analyzed for the solutions to the elliptic counterpart of problem (1.1) for general 
structures of p-Laplacian type, again taking into account the original potential 
approach in [32] 2 . 

In this paper, we will present a series of up-to-the-boundary regularity results 
aimed at giving, on one hand, an extension of the classical Lebesgue theory devel- 
oped in [8] to a very general class of rearrangement invariant function spaces. For 
instance we will cover the case of Lorentz spaces, that do not seem to be treatable 
with the currently available methods. On the other hand, we will give a set of 
techniques allowing for treating domains with very poorly regular boundaries. In 
this respect we also extend to the parabolic case the recent results obtained in [40]. 

Indeed, since we are dealing with precise estimates up to the boundary, it is 
necessary to take some regularity condition for = ^ x (— T, 0). For this, we 
assume that the domain has only to satisfy a capacity density condition. Namely, 
in order to achieve the main results of the present paper, the complement of the 
set £1 is required to be a thick domain (see Section 2.4 below); and such assumption 
is both very natural (as, for instance, in the classic general Sobolev embedding 
in capacity density settings; and, also, in the elliptic framework studied in [40]), 
and very mild (it is fulfilled by any domain with Lipschitz boundary, or even of 

^We refer to Section 2.2 for the definitions of the involved function spaces in the parabolic 
setting. 

2 Although - as stated in [40], too - the possibility of extending such local results to global ones 
was already mentioned without proof in [32] and [34] . 




locally in Qt 
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corkscrew type; also, in K 2 by any proper simply connected domain). As a further 
matter, this condition is essentially sharp for higher integrability results, as pointed 
out by Kilpelainen and Koskela in Remark 3.3 in the paper [23], where the authors 
developed a first analysis in global Lebesgue regularity of the gradient of solutions 
to degenerate elliptic equations in general thick domains. 

All in all, we extend the potential theoretic strategy recently developed by Min- 
gione in [32, 33] to the parabolic framework up to the boundary, combining it with 
global higher integrability results (in Theorems 3.5 in Section 3), to obtain general 
gradient estimates on level sets (see Section 4), which in turn imply the following 

Theorem 1.2. Let Q be such that H n \ is a thick domain with constants cq and 
po and assume (1.2). Let f € £(7, q)(p,x) with 



where c is a constant depending only on n, q, L, u, 7, cq, po> diam(fi). 

Of course, by choosing q = 00, we can also deduce regularity results on the 
Marcinkiewicz scale. 

Remark 1.3. After stating Theorem 1.2, we want to make a more precise comparison 
between our results and the basic results in [8], where the authors treat the case of 
Lebesgue spaces; i.e., they take q = 7 in the statement above thereby obtaining the 
corresponding conclusion. The results in Theorem 1.2 is not a direct extension of 
the ones in [8], as in this paper no regularity assumptions are considered on £1. On 
the other hand, the results in the present paper are global but they are obtained 
by covering arguments via estimates that are local in nature as it is easy to see 
that, using the methods developed here, one can obtain local estimates up to the 
boundary and not only in the interior. This corresponds to replace, in Theorem 
1.2, fir by cylinders of the type (B(xo, R) D fi) x (-T, 0), where xq € dQ. This is a 
valued consequence of the techniques approached here and the necessity of assuming 
the considered boundary capacity conditions stems exactly from this aspect of the 
involved techniques. Now, whereas it is likely that some Lorentz space estimates 
can be also obtained under milder assumption on Q, it is worth mentioning that 
the methods here can be also used to get global estimates in Morrey and Lorentz- 
Morrey spaces (as those obtained in the interior case in [14, 32]). 

The proof of Theorem 1.2 relies on the fact that the integrability of the spatial 
gradient of the solutions to problem (1.1) is linked to a suitable choice of potential 



Then the solution u E L 





Moreover, the following estimate holds 
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operators. In this sense, one of the key-points will consist into obtaining a decay 
estimate that involves the level sets of the parabolic Hardy-Littlewood maximal 
operator Mq of \Du\ in term of those of another maximal operator M* of the 
datum /, up to a correction term which is negligible when considering the higher 
regularity. Indeed, we will prove an estimate of the following type 

(1.5) \{Mq(\Du\) > 5A}| < J-|{M *(|LH) > A}| + |{[M*(|/|)] > A}|, 

for every A suitably large, and in which S » 1 is a constant to be chosen, and the 
exponent %i > 1 is related to the higher integrability theory. 

The estimate in (1.5), whose precise version is given by forthcoming formula 
(4.12), is fairly general and it will be relevant to deducing the Lorentz estimates 
stated in our theorems, also including the borderline case. In order to obtain such 
level sets estimate, we apply a parabolic covering lemma of C alder on- Zygmund type 
(see Lemma 2.12 below) together with the higher integrability theory which we will 
extend up to the boundary (see Theorem 3.5 below). Therefore, we will work on 
basic estimates of the solutions u to (1.1) in comparison to the solutions w to 
the corresponding homogeneous problem (see Section 3), by exploiting very recent 
contributions in the parabolic framework given in [27]. Clearly, the situation is 
complicated by the necessity to work in a global setting. 

Analogously, we will show how to establish a modified version of the level sets 
estimate in (1.5) by replacing Du by u there. This will permit us to recover higher 
regularity estimates in the Lorentz scale for the solutions u themselves, as stated 
in the following 

Theorem 1.4. Let fl be such that W b \ O is a thick domain with constants Co and 
po and let f G L(j, q)(£lT) with 

N 

1 < 7 < — and < q < oo. 
Then the solution u € L l {— T, 0; Wq to (1.1) satisfies 

ueL [N=*j' 9 ) m ^ T - 

Moreover, the following estimate holds 

IMI^i,)^) ^ c ll/IU(7,*)(nr). 

where c is a constant depending only on n, q, L, v, 7, cq, po, diam(O). 

It is worth pointing out that, although in the elliptic case the regularity of u 
can be recovered by plainly combining the regularity of the gradient Du with the 
classic Sobolev embeddings, here we need to work in a separate way, by means 
of sharp estimates also involving different maximal operators in (1.5), and higher 
regularity results up to the boundary for the solutions u. 



Finally, we recall that all the results we obtained in the present paper hold for the 
very weak solutions given by the approximation method described in forthcoming 
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Section 2.1. It would be interesting to understand whether these results can be 
extended to some other notion of solutions. In this respect, a positive answer can 
be given when dealing with notions of solutions to measure data problems holding 
uniqueness in case of integrable data, as in case of the renormalized solutions in [38] 
and [39]; see [7] for a first definition of renormalized solutions in this framework, 
and also [16]. 

Remark 1.5. In the present paper we are dealing with the case of operators with 
linear growth while in other papers, as for instance [32], level sets techniques are 
used to approach similar regularity results for p-Laplacian type operators, that is 
for equations of the type — div (\Du\ p ~ 2 Du) = f. When considering the parabolic 
analog 

u t -A\-v{\Du\ p ~ 2 Du) = /, 

the approach of the present paper is not sufficient if used alone, and has to be 
modified and extended using the concept of intrinsic geometry first pioneered by 
DiBenedetto to overcome the lack of homogeneity of the evolutionary p-Laplacian 
operator when p / 2 (see for instance [1] for a discussion of this concept in the 
relevant context). This is the object of forthcoming investigations ([6]). 

The paper is organized as follows. 

In Section 2, we fix notation; we give some details on the structure of the problem 
and we briefly recall the definitions and a few basic properties of the spaces and 
the operators we deal with, also discussing the capacity density condition and some 
properties of Sobolev functions in p-thick domains; we recall a parabolic covering 
lemma of C alder on- Zygmund type. 

In Section 3, we state and prove global higher regularity results for the solutions 
to the homogeneous analog of problem (1.1), and other preliminary results. 

In Section 4, we prove the results stated in Theorem 1.2, and the Lorentz space 
estimates for the solutions u. 

2. Preliminaries 

Throughout the paper we follow the usual convention of denoting by c a general 
positive constant always greater or equal than one that may vary from line to line. 
Relevant dependencies will be emphasized by using parentheses; special constants 
will be denoted by cq, c±, C2,.... 

As usual, we denote by 

B R (x ) = B(x ; R) := {x G R n : \x - x \ < R} 

the open ball centered in xq 6 M. n with radius R > 0. When not important 
and clear from the context, we shall use the shorter notation Br = Br{xo). Let 
(x ,t ) € R n+1 , we denote by 

(2.1) Qr(x , t ) = Q{x , t ; R) := B{x ; R) x (t - R 2 ,t ) 

the open parabolic cylinder with height R 2 and having a ball Br as horizontal slice. 
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Finally, we mean <9 par ^T = Aat^T U Qo, where <9i a t^T := dCl x (— T, 0) is the 
lateral boundary and Qq = Q x {— T} is the initial boundary. 

2.1. Solvability of the problem. A measurable function u is a very weak solution 
to (1.1) if ue L^-^OjWo 1 ' 1 ^)) and 

(2.2) - / u(f> t dxdt+ / (a(x,t,Du),D(f>)dxdt = / f(f)dxdt, 

Jn T Jn T Jn T 

for any <j) £ Cq°(Qt)- Also, while the lateral boundary condition can be formu- 
lated by prescribing the belonging of u to L 1 (— T, 0; W ' (fl)), the initial boundary 
condition is understood in the L 1 -sense, that is 

1 r-T+h r 

(2.3) -/ / \u(x,t)\dxdt -> as/i->0. 

As usual in the parabolic setting, one can provide a convenient "slice-wise" re- 
formulation of equality (2.2) by mean of Steklov average. Indeed, for h > and 
t € [— T, 0), we can define 

( 1 rt+h 



u h (x,t) :-- 



I rt+n 

- / u(x, i)dt iit + h<0 
if t + h > 



and the following equality 

(2.4) / [(u h ) t <f>+{[a(x,t,Du)] h ,D<j>)]dx= I f h <f>6x 



holds for any <j> G Wo' 2 (n) and for a. e. i € (— T, 0). This formulation allows to use 
the solution u, which enjoys only very weak regularity properties in time, as a test 
function in (2.4), and therefore throughout all the paper we shall proceed formally 
in this sense, referring for instance to [4] for the details in similar situations. 

The existence of such solution is obtained using an approximation method; this 
means that we shall consider a special class of solutions called SOLA (Solutions 
Obtained by Limits of Approximations). In particular, one considers a sequence 
of bounded function {fk} C L°°(Q,t) such that fk — ¥ f in L 1 (^t) as k — > oo. 
Then, by standard monotonicity arguments, for each k fixed, there exists a unique 
solution 

u k € C7°([-T, 0];L 2 (n)) n L 2 (-T, 0; ^(ft)) 

and 

to the Cauchy-Dirichlet problem 



(2.5) 



(u k ) t - div a(x, t, Du k ) = f k in tt T 

U k = on dpar^T- 



Finally, thanks to the arguments in [9, 8], we can pass to the limit in the problem 
above and to prove the existence of a function u € L l (—T, 0; Wq' 1 ^)) which solves 
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(1.1) in the distributional sense. For the rest of this paper we understand by {uk} 
the sequence obtained by solving (2.5) with 

(2.6) fk(x,t) :=max{ -k, min{/(x,i), k}}, feel. 

Notice that, as stated in Remark 1.1, this choice is not restrictive, since in this class 
of solutions uniqueness holds, in the sense that considering another approximating 
sequence {fk} would lead to the same limit u. 

2.2. Relevant parabolic function spaces. Let £It be the space time cylinder 
f2 x (— T, 0), with being a bounded open set in R n , n > 2 and T > 0. 

Fix q G (0, oo). A measurable map g : Q,t R fc belongs to the Lorentz space 
L(-y, g , )(flr) with 1 < 7 < 00 if and only if 

f°° 1 riA 

( 2 - 7 ) bWlfaxnT) := q J, G Ut : \9(x,t)\ > A}|)^ — < +00. 

In the case q = 00, the Lorentz space £(7,00) with 7 G [1, 00) is the so-called 
Marcinkiewicz space and it is usually denoted by M^^t)- A measurable map 
g : Q T — > ]R fe belongs to M^^t) if an d only if 

(2-8) llsllXnpi,,) = ll^lli( 7 ,oo)(n T ) := supA 7 |{(x,t) g n T : \g(x,t)\ > X}\ < +00. 

Note that by Fubini's Theorem one can see that, for any 7 G [1, 00), £(7, 7) (fir) 
is nothing than the classic Lebesgue space L~ ( (£1t)- Also, by Fatou's Lemma, one 
can see that the functionals || • ||l( 7)IJ )(o, t ) defined in (2.7)-(2.8) are lower semi- 
continuous with respect to the a. e. convergence; see [5, Section 3.2]. 

Moreover, the spaces defined above enjoy Holder type inequalities. In particular, 
below we state a standard inequality for the Marcinkiewicz spaces A1 7 (0) in the 
form we will need it in the following of the paper. 

Lemma 2.1. Let E C R n+1 be a measurable set and let g G J^V{E) with 7 > 1. 
Then, for any q G [1,7), g G L q (E) and 

( 7 \ ^ 

\\9\\^{E) < [^^) \ E \ q ~<\\9\\m^e)- 

For details and results about the theory of Lorentz spaces, we refer the interested 
reader to [2, 20]. 

2.3. Parabolic maximal operators. For any measurable function g, the re- 
stricted fractional maximal operator M* g , with a G [0, iV], relative to a parabolic 
cylinder Qo, is defined by 

M* iQo (g)(x,t) := sup \Q\%1 \g{y,T)\dydT, 

QCQo, {x,t)eQ JQ 

where the cylinders Q have sides parallel to those of Qo- 
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The boundedness of maximal operators in Marcinkiewicz spaces is classical (see, 
for instance, [20]); i.e., for any g 6 L^(Qq), 

(2.9) \{(x,t) € Q : Ml Qo {g){x,t) > X}\ < [ \g\^6xdt 

holds for every A > and 7 > 1; the constant c depending only on n and 7. 

More in general, a standard embedding result for the maximal functions in 
Lorentz spaces does hold, as given by the following theorem, whose proof follows 
by the Marcinkiewicz theorem together with standard sublinear interpolation. 

Theorem 2.2. ([32, Theorem 7]). Let a in [0, N), and 7 > 1 be such that cry < N. 
Let Q C be a parabolic cylinder and q 6 (0, 00]. Then for every measurable 

function g in L(-y,q)(Q) it holds 

\\ M a,Q(9)\\ L r_Nj_ jq y Q) < c lbllL( 7 ,g)(Q), 

where c is a constant depending only on a, 7, n and q. 

Also, in the borderline case 7 = 1, we have the following 

Theorem 2.3. ([5, Lemma 4.11 and Theorem 4.12]). Let a € [0,N), Q C R n+1 
be a parabolic cylinder and consider the concentric parabolic cylinder o~Q scaled 
by a factor a > 1. Then there exists a constant c = c(n,a,o~) such that, for any 
measurable function g in o~Q, it holds 

H M ^)II L7V ^ (Q) ^ c l Q l 1-f ll^l5(.a)ll^l^oi( S) - 

2.4. Capacity density condition and thick domains. In the following we re- 
call the definition and some properties of the p-thick domains. We need first to 
recall the definition of p-capacity of sets. 

Definition 2.4. Let if C be a compact set. For any p > 1 the variational 
p- capacity of K on Q is defined by 

cap p (K,0) := inf I J \Dip\ p dx : <p G C™{£l),y > 1 K \ . 

As customary, the definition above can be extended to generic open sets ACS] 
by taking the supremum over the capacities of the compact sets contained in A. In 
the case of balls, the definition simplifies considerably and we have the following 
explicit formula for B p on B<i p 

civp p (B p ~,B 2p ) = cp n ~ p , 

where c is a positive constant depending only on n and p. For further details about 
the variational p-capacity, we refer for instance to Chapter 2 in [21]. 

Definition 2.5. A set E C W 1 is uniformly p-thick (equivalently, satisfies the p- 
capacity uniform thickness condition) if there exist constants cq and po > such 
that 

(2.10) cap p (Er\B p (x),B 2p (x)) > c c&p p (B p (x), B 2p (x)) = c(n,p, c ) p n ~ p , 
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for any x G E and for any < p < po- 

For the sake of simplicity, from now on we shall use the short notation thick to 
denote a 2-thick domain. 

As stated in the introduction, one can prove that domains satisfying Defini- 
tion 2.5 include those with Lipschitz boundaries or even those that satisfy a uniform 
corkscrew condition. Furthermore, the estimate in (2.10) remains valid for balls 
centered outside a uniformly p-thick domain (and near the boundary), as stated in 
the next lemma; see [37, Lemma 3.8]. This will be useful since we are going to deal 
with bounded sets whose complementary is thick. 

Lemma 2.6. Let Q be a bounded domain such that W n \ fi is uniformly p-thick 
with constants cq and po, and let y G Q be such that B p /§(y) \ / 1 Then there 
exists a constant c\ = ci(n,p, pq,cq) such that 

c&Y> p (B p/A (y) \ Q, B p/2 (y)) > ci c&p p (B p/4 (y),B p/2 {y)) = cp n ~ p . 

Also, a uniform p-thick domain enjoys a self-improving property, as stated in 
the following theorem due to Lewis [28] (see also Section 8 of [30] for a complete 
survey on the boundary regularity). 

Theorem 2.7. Let Q be a bounded domain such that W 1 \ Cl is uniformly p-thick 
with constants cq and po, and 1 < p < n. Then there exists q = q(n,p,co) G 
(l,p) such that R n \Q is uniformly q-thick with constant c 2 = c 2 (n,p,co) and pq. 
Moreover, q can be chosen near p so that q G (np/(n +p),p). Thus, we have 

cap q (B p/i (y) \n,B p/2 (y)) > c 2 c&p q (B p/i (y), B p/2 (y)) > c{n,p,c )p n ~ q 

for y G £1 be such that B p / 6 (y) \ Q ^ 0. 

Remark 2.8. On the other hand a p-thick domain is r-thick for every r > p. Indeed, 
for x e E, p < p and ip £ C^(B 2p (x)), <p > ^ En B t^j , by the Holder inequality, 
we have 

f \Dp\ p dx < c{n,p,r) ( [ \D^\ r dxY p 71 ^) . 
Jb 2p \Jb 2p ) 

Taking the infimum for such ip we get 

c(n,p, c ) p n ~ p < cap p (£ n B~ p , B 2p ) < cp n ( 1 "?) [cap r (£ n B~ p , B 2p )] ? 
and r-thickness follows by taking the r/p-power. 

We conclude this section by recalling some capacity Sobolev type inequalities; 
see for example [23, Lemma 3.3] or [30, Lemma 8.11 and Remark 8.14]. 

Lemma 2.9. Let w G W l ' p (B 2p ) and K p (w) = {x G ~B~ p : w(x) = 0}. Define 
k = n/{n — p) if p < n and k = 2 if p > n. Then there exists a positive constant 
c = c(n,p) such that 

| I \iu\ Kp dx | < ( - f \Dw\ p dx | 

\Jb 2p J ~ \cap p (K p (w),B 2p ) J B2p J 
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The previous inequality also holds in the case stated below, with p-powers on 
both sides; that is, we have the following Poincare type inequality 

Lemma 2.10. Let w G W 1 ' p (B 2p ) and K p (w) as before. Then there exists a 
positive constant c = c(n,p) such that 

(2.11) / M p dx< C f \Dw\ p dx. 

JB 2p cap p {K p (w),B 2p ) J B2p 

Now a parabolic Sobolev inequality; its proof requires only slight changes with 
respect to the proof of Lemma 3.17 in [37]. Essentially, it is enough to apply 
Sobolev's inequality Lemma 2.9 slicewise to the continuation to zero of w in Q p / 2 \ 
VLt- Note moreover that the exponent q appearing therein is the same present 
in Theorem 2.7, for p = 2, and that here it is unavoidable such a self-improving 
property. 

Lemma 2.11. Let 0, C R n be a bounded domain such that R n \Q is uniformly 
thick with constants Co and po, and let Q p /2 '■= B p /2( x i) x (*i — (/>/2) 2 , ti) be a 
parabolic cylinder such that B p / e (xi) \ H / 0. Let moreover w be a function such 
that 

w £ L 2 {Q p/2 nn T ), DweL 2 (Q p/2 nn T ), 

w eL°°(t 1 -(p/2) 2 ,t 1 ;L 2 (B p/2 (x 1 ))) 

and w = on Q p / 2 H Siat^T- Then there exist an exponent q = q(n,co) < 2 and a 
constant c = c(n,co, po) such that 



(2.12) / \w\ 2 dxdt<c — — - \Dw\ q dxdt 



Qp/2 1 JQ p/2 nn T VI Qp/2 1 JQ p/2 nn 



x ( sup / |u?(-, r)| 2 dx 

vrG(ti-(p/2) 2 ,ti)n(-T,o) JB p/2 nn 



2 
n + 2 



and 



(2.13) / \w\ 2{1+ ^ dxdt <c I \Dw\ 2 dxdt 

JQ P / 2 nn T JQ p/2 nn T 



x I sup / \w(-,t)\ 2 dx 

\re(ii-(p/2) 2 ,ti)n(-T,o) JB p/2 nn 

2.5. A parabolic covering lemma. Below we present a technical lemma (see 
forthcoming Lemma 2.12) which is nothing that a special version of the classical 
Calderon-Zygmund-Krylov-Safonov decomposition. It will allow us to work with 
parabolic cylinders that have as horizontal slice a ball instead of a cube. We would 
recall that a parabolic version of the Calderon-Zygmund-Krylov-Safonov covering 
lemma was proved by one of the authors in [5] , following the elliptic analog (which 
could be found, for instance, in [11, Lemma 1.2]). Precisely, in [5, Proposition 2.1] 
it is shown that a covering of a set by parabolic cylinders can be "arranged" in 
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an optimal disjointed way so that in each of the relative horizontal slice cube the 
density is considerably small. 

Since in the present paper we are interesting into analyzing some regularity 
properties of problem (1.1) up to the boundary, we need to work with parabolic 
cylinders having a ball as horizontal slice. For this, in order to prove Lemma 2.12 
below, we will extend its elliptic counterpart given by Lemma 3.1 in [40] to a 
parabolic framework. 

Lemma 2.12. Let Q R C R n+1 . Assume that X C y C Q R are measurable sets 
such that there exists < 5 < 1 for which the following properties hold: 

(i) |*|<<f|Qfl|; 

(ii) for all (x, t) £ Q R and for all p€(0,R] if \X PI Q p (x, t) | > 5 \ Q p (x, t) \ then 

Q P (x,t)nQ R cy. 

Then we have 



where C3 = 03(71). 

Proof. Since (i) holds, for almost every (x, t) € X there exists a p x t < R so 
that \X n Q Px . t (x,t)\ = 5\Q Px J and \X n Q P (x,t)\ < S\Q p \ for all p x<t < p < R. 
Then, by the Vitali covering lemma we can extract a sequence of disjoint cylinders 



{Qpx k ,t k ( x k,tk)} k such that X C (\J k Q5p Xki t k (x k ,t k )j n Q R . In view of the choice 



X n Q 5pxkttk (x k ,t k )\ < 6\Q 5pxk<t J = 5 n+2 5\Q Pxk J = 5 n+2 \X n Q Pxk , tk (x k ,t k )\. 



x\ <c 3 6\y\ 



of Q 



(x k ,t k ), we infer that 



Moreover, there exists a constant c{n) such that 



s, 



(x k ,t k )\ < c(n)\Q t 



(xk,t k ) n Q 



since (x k ,t k ) G and /j Xfejffc < R. 
All in all, we have 




k 



< 



k 



< 



c(n)5 n+2 Y,S\Q P , k , tk (x k ,t k )nQ 



R 



k 



< 



c(n)5 n+2 5\y\, 



where we used (ii) in the last inequality. This concludes the proof. 



□ 
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3. (Higher) regularity for the homogeneous problem and comparison 

RESULTS 

One of the key-points in the proofs of the main results in the present paper relies 
on obtaining both local interior and boundary comparison estimates. 

Let us first handle the interior ones by considering the following homogeneous 
Cauchy-Dirichlet problem. 

{wt — div a(x, t, Dw) = in Qr 
W = U On dparQfi, 

with Qr = Qfi(a;o,io) C and u being the unique solution to the regularized 
problem (2.5). We recall some basic results from the higher integrability theory of 
Gehring. For the proofs we refer to [29]; see also [5, Theorem 5.5]. 

Theorem 3.1. Letw G C°([t -R 2 , to] I L 2 (B R (x )))nL 2 (t -R 2 , t ; W lfl (B R (x ))) 
be a weak solution to the parabolic equation (3.1)i. Then there exists xi = Xi( n ) L, v) 
> 1 such that Dw G L 2 ^(Qr) and for any q G (0, 2] it holds 



(3.1) 



i 

(3.2) (/ \Dw\ 2xi dxdt J <c(f |Diu| ff da;dt ) 

\jQ R /2 J \JQr ) 

Moreover, for any q G (0, 2] and for any x > 1 it holds 



(3.3) (/ \wf x Axdt\ <c(-f \w\ q dxdt 

The constants c in (3.2) and (3.3) depend only on n,L,u and q. 

Concerning local interior comparison estimates we can use those established in 
[5, Lemma 5.2] or [4, Lemma 6.4]. 

Lemma 3.2. Let u G C°([-T, 0]; L 2 (tt)) n L 2 (-T, 0; W 1,2 (fi)) be the unique solu- 
tion to the regularized problem (2.5). Moreover, let w be the unique solution of the 
Cauchy-Dirichlet problem (3.1). Then there exists a constant c = c{n,u) such that 

/ (R^\u-w\ + \Du- Dw\)dxdt < cR |/|dx<tt. 
JQr JQr 

For the local boundary estimates we consider Qr = Qr(xo, to) C M. n+1 such that 
it intersects the lateral boundary of Qy, that is 

Qr n (dn x (-t, o)) = Qr n d ut n T ? 0. 

This will be enough for our purposes, see Remark 4.1 at the beginning of Sec- 
tion 4. Then, we consider the unique solution to the following homogeneous 
Cauchy-Dirichlet problem 

| w t - div a(x, t,Dw) = in Q R n fir 

I w = u on d par (Q R n Or). 
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Here we need a higher integrability result for w, the counterpart of Theorem 3.1 
up to the boundary, relying on some reverse Holder inequalities. This kind of 
boundary higher integrability is already present in the literature (see for instance 
[37, Theorem 4.7], and also [3, 10, 36] for more regular boundaries), but we prefer 
to propose an independent statement (and subsequently the relative proof) fitting 
our context. 

For this, we need a couple of technical lemmata, the first of which is an appro- 
priate version of the Gehring lemma, see [18, Proposition 1.3]: 

Lemma 3.3. Let A C R n+1 an open set and let g : A — >• M, k be an integrable map 
such that 

2 

f \g\ 2 dx dt < c ( J- \g\ q dxdt ) * +■&+ \g\ 2 dxdt, 
Jq p / 6 \JQp / JQp 

for some constant q < 2 and for all Q p C A. Then there exists a constant $o = 
•&o(n,q) such that if $ < $0 then g E L^(A) with xi > 1 depending on n,c,q,-d; 
moreover 



(/ \g\ 2xi dxdt ) 1 < cl \g\ 2 dxdt 



for all Q p C A. The constant c also depends on n,c, q,$. 

The second one encodes the self-improving character of reverse Holder inequali- 
ties. The proof follows [19, Remark 6.12] and uses a covering argument. 

Lemma 3.4. Let A, g as in Lemma 3.3 and let 



(/ \g\ 2xi dxdt) X1 < cl \g\ 2 dxdt 

\JQn/2 / JQp 



ip/2 

hold whenever Q p C A is a cylinder, where xi > 1 and c > 1. Then, for every 
a € (0,2], there exists a constant c = c(n,o~,c) such that 

(J \g\ 2xi dxdt^ 2 * 1 <c(j \g\ a dxdty 

for every Q p C A. 

Now, we are in a position to prove the following 

Theorem 3.5. Let be a bounded domain such that W 1 \ Q satisfies a uniform 
thick condition with constants cq and po, and let w be the unique solution to (3.4). 
Then there exists a constants xi = Xi( n > ^ L, po, cq) > 1 such that for every q > 
there holds 

1 

( 3 - 5 ) I m 1 I \Dw\ 2xi dxdt) <c(-^—[ \Dw\ q dxdt). 

\\Qr/2\ J Q R/2 nn T J \\QR\JQ R nn T J 
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for a constant c = c(n, u, L, po,co, q). Moreover we also have 

( 3 - 6 ) I f \w\ 2x dxdt ] < c(—^— [ \w\ q dxdt 

\\Qr/2\ JQ R/2 nn T J \\QR\jQ R nn T 

for any x = x( n i v i L, Pq, cq) > 1 and the same constant as above. 

Proof. The goal here is proving that if we consider the continuation of w to zero 
(denoted by w) in Qr\£It, then for every parabolic cylinder Q p = Q p (xi,ti) C Qr 
we have a reverse Holder's inequality 

(3.7) (J \Dw\ 2 dxdty < c(J \Dw\ q dxdty + ^(j |ZM;| 2 dxd^ 2 , 

for some q < 2 and for a constant depending on n, u, L, po, cq, #o> where i?o is the 
constant of Lemma 3.3 corresponding to q. Therefore we could apply the Gehring 
lemma in the form of Lemma 3.3, and then Lemma 3.4, to infer (3.5), after again 
restricting w to Qr n Qt- The point here is that we are not interested in the 
boundary higher integrability of w in the part of the same boundary inside CIt, but 
just in the higher integrability on cylinders centered in d\ at ^T, where w = u = 0. 
In a very rough sense, this is similar to interior higher integrability. Hence now we 
fix a parabolic cylinder Q p C Qr; note that due to the assumption Qr n f^o = 0, 
we also have Q p n f^o = 0- 

The scheme of the proof is the following: first we consider the case where -Bp/6 \ 
/ fl. Here we first use a Caccioppoli type inequality to, roughly speaking, 
control the gradient Dw over -B p/ / 4 n £It with the integral of w over B p / 2 fl £It- 
By applying the Sobolev inequality over B p j2 fl Ot, we now need to estimate the 
supremum in (3.14) in terms of \Dw\ and this is again achieved with the help of a 
Caccioppoli type inequality and the Poincare capacity inequality in Lemma 2.10. 
Finally, we get (3.7) by narrowing the set on the left-hand side. On the other hand, 
if -Bp/6 \ O = 0, estimate (3.7) between -B p / 8 and -B p / 6 is just an interior one, and 
then we can enlarge the domain on the right-hand side. 

Suppose now hence B p /e \ 7^ (but clearly also -Bp/6 H ^ ^ 0, otherwise 
there would be nothing to prove) and consider a cut-off function <p € Cq°(Qr), 
< ip < 1 such that (p = 1 on Q p u, suppy? C Q p /2 an d \Dip\ + p\ipt\ < c/p. Note 
that from now on all the cylinders with radii multiple of p will have vertex (xi,t\). 
Up to a standard regularization of w in time by Steklov averages, we use as a test 
function in (3.4) the map (j> := wtp £ £ 1q t , where, once fixed t 6 (h - (p/2) 2 ,ti), 
Q = Ce,r(t) G M /1 '°°(IR) is a piecewise linear continuous map, for e < (t± — r)/2, 
such that ^ e = 1 in (— 00, r] and Ce = in [r + e, +00). Note that for a. e. 
t e (tx - (p/2) 2 ,t 1 ), we have cp(-,t) € Wq 2 (Q) and thus can be used in (2.4). We 
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therefore have, integrating over (t\ — (p/2) 2 ,t±), 

(3.8) / wtW(p 2 ( £ dx dt + / (a(x,t,Dw),Dw)tf 2 ( £ dxdt 

Jq p /2^t JQ p /2nn T 



-2 / (a(x,t, Dw), D(p)wip( £ dx dt. 



Now, we treat the first term in the following way, using integration by parts and 
recalling that ip( £ = on (B p / 2 Pi f2) x {ti — (p/2) 2 ,t{\, 



w t w<p 2 ( £ dxdt = ^ / (\w\ 2 ) t ip 2 ( £ dxdt 
p/2 nn T 2 JQ p / 2 nn T 



\w\ 2 (ip 2 ( £ ) t dx dt 
Q P / 2 nn T 



(3.9) 



f ,2 1 f T+£ f ,2 2 

/ (ptftCe dx dt + - + / </? dxdt. 



2 p/2 nQ T z A iB p/2 r 

Letting e \ and then taking the (essential) supremum with respect to r € 
(ti — (p/2) 2 ,ti) we therefore get 

(3.10) sup [ wtwip 2 dx dt > — c f J — j-dxdt 

re(ti-(p/2) 2 ,ti) JQ p /2,r^T JQ P /2nn T P 

+ - sup / [\w\ 2 ip 2 ]{-,T)dx. 
1 re(ti-(/9/2) 2 ,ti)iB p/2 nn 

where we denoted for shortness with Q P /2, T the set Q p /2 H l.r i<r }. The second term 
in the left-hand side of (3.8) is estimated from below 

(3.11) / (a(x,t,Dw),Dw)ip 2 ( £ dxdt > - \Dw\ 2 (p 2 dxdt 
jQ p/2 nn T 2 ^s p/2 , T nn T 

with c = c{u, L), having at hand the following monotonicity condition, that can be 
deduced by the assumptions in (1.2): 

(a(x,t, z), z) > 2^ z \ 2 ^ or an y ( x ify ^ and any z € Il n . 

For the latter term in the right-hand side of (3.8) we have, using again the 
growth condition (1.2)2 together with the Young inequality, 



(a(x, t, Dw), Dip)wipCe dx dt 



(3.12) <~ f \Dw\ 2 ip 2 dxdt+ C ^' 2 L ^ f \w\ 2 dxdt. 

4 Jn T nQ p/2 ,r+e P JQ P / 2 r\n T 

To conclude we put (3.9), (3.11) and (3.12) into (3.8), we let e \ 0, we perform 
some algebraic manipulations and we take the supremum with respect to r € 
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(fi - (p/2)Vi). Using also (3.10) finally we get 

sup / \w(-,r)\ 2 dx + / \Dw\ 2 dxdt 

re(ti-(p/4) 2 ,ii) JB p/4 nn JQ p/i c\Q, T 

(3.13) <£y^) I' \ w \ 2 dxdt, 

P JQ p/2 nn T 

recalling that (p = 1 on Qp/4. Note that also a completely similar estimate, where 
the cylinder Qp/2 appears on the left-hand side, and the cylinder Q p on the right- 
hand side, follows by straightforward changes in the proof above. 

We recall now we are supposing B p /§ \ f2 ^ 0; therefore we can use the Sobolev 
inequality in Lemma 2.11 in the right-hand side of (3.13) and, after dividing by 
| Qp/2 1, we infer 

\Dw?dx dt 



|2p/4l JQ p/4 nn T 

2n 

/if , \^+2) 

(3.14) < c I — -/ \Dw\ q dxdt) 

VI Qp/2 1 JQ p/2 nn T J 

. 2 

sup ttt — i / K-' T )I dx 

re{t 1 -( P /2) 2 ,t 1 ) \^p/2\ JB p/2 nn J 

where c depends upon n, v, L, po> c o- Finally, we estimate the supremum on the 
right-hand side in the following way: for a. e. r € (ii — (p/2) 2 ,ii) using (3.13) in 
its version over Q p /2 and Q p and Poincare's inequality (2.11) slice-wise, for p = 2, 
we have 

\w(-,t)\ 2 dx < — \w\ 2 dxdt 
B p/2 nn P JQ p nn T 

4 f 1 7— — — -r / |Du;| 2 dxdt 



p 2 J tl _p2 cap 2 (ir p/2 (u;),Sp) ,/ B 



where K p / 2 {w) = {x € -Bp/2 : iw(x) = 0}. We recall that we extended w to zero in 
Qr \ Qt', note this gives w £ W 1,2 {Br) for a. e. i since w = on di at Q,T- Now, 
since K p / 2 (w) D B p / 2 \ Q and -Bp/ 6 \ £1 ^ we can apply Lemma 2.6 to get 

cap 2 (K p/2 (w),B p ) > cicwp 2 {B p/2 ,B p ) = cp n ~ 2 , 

and therefore, again dividing by |Q P /2|> 



SU P Lo — i / |^(^ r )| 2 dx < — r / \Dw\ 2 dxdt. 

re(ti-(p/2) 2 ,ti) I Qp/2 1 JB p/2 nn I Qp/2 1 JQ p nn T 

Putting the latter estimate into (3.14) and using the Young inequality with conju- 
gate exponents (n + 2)/2 and (n + 2)/n, with "#o being the constant appearing in 
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Lemma 3.3 corresponding to the choice of q following from Lemma 2.11, we infer 
\Dw\ 2 dx dt 



Qp/i\ JQ p/i nn T 



<c|— — - \Dw\ g dxdt) — — ,/ \Dw\ 2 dxdt 

I Sp/2 1 JQ p/2 nn T J V I Qp/2 1 JQ p nn T 

2 
9 



- C ^°) ( Tn — i / \Dw\ q dxdt 
\Qp/2\ JQ p/2 nn T 



■(—') 



2 1 



+ ( "77777 ) — T / I^Pdxdt. 



viooy |Q P |7 Sp n^ T 
The inequality above plainly yields (3.7) in the case B p / G \ ^ 
The case B p /6 \ ^2 = is easier, since now an estimate 



J- \Dw\ 2 dxdt Y < c ( f \Dw\ q dxdt) 



q possibly different from the one appearing in (3.14), is an interior one and can 
be deduced following [18] or [24]. At this point we just enlarge the integral on the 
right-hand side to get (3.7). Note that actually higher integrability estimate (3.5) 
holds for every couple of cylinders Q p /4, Q p /2 such that Q p C Qr, and we shall 
use this fact in few lines. 

Finally, the estimate (3.6) for w is a straightforward consequence of boundary 
sup-estimate for parabolic problems with null boundary datum. For this, we refer 
for instance to [29, Theorem 6.30] in a slightly different, linear setting. For the 
sake of the reader we also propose in the following lemma a proof adapted to our 
framework. □ 

Lemma 3.6. Let w be the solution of the Cauchy-Dirichlet problem (3.4). Then 
for every q > it holds 

(3.15) sup \w\ <c(— - — | I \w\ q dxdt 

Q R/2 rin T \\Qr\ JQ R nn T 



for a constant c depending n, u, L, cq, po, q. 

Proof. Consider, for R/2 < p < R and ft£l, the quantities 

U(h,p):= / \(w — h) + \ 2 dx dt, 



(3.16) 

v(h, P ) -.= \Q p nn T n{w > h}\ =: \A(h,p)\; 
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here s+ is the positive parts of s; here all the cylinders share the same center (xo, to)- 
For p\ < p2 both belonging to [R/2,R] and with h < k, we have the relations 

V(k, Pl ) < —±— U(h,p 2 ), 

(3.17) U(k, Pl )<c- 1 U(h, P2 )[V{h,p 2 )}^. 

[P2 ~ Pi) 2 

The first inequality plainly follows by the involved definitions. Indeed, 

(3.18) V(k, Pl )(k - hf < / (w-h)\dxdt< (w - h)\dxdt. 

J A(k,pi) J A(h, Pl ) 

For the second one, we need some work. First of all, minor modifications in the 
proof of Caccioppoli's inequality (3.13) give for all i € IN 



sup 

re(<o-(pi+P2) 2 /4,to) JB iPl+P2)/2 nn 



/ \(w - k) + (-,T)\ 2 dx 

JB, n , _i_„„-, / 9 nf2 



2 



+ / \D(w - k) + \ dxdt 

(3.19) < - C ^ ,L \ 2 f \{w -k) + \ 2 dxdt 

{P2-pi) JQ P2 nn T 

for every h € R; in particular we need to choose here the test function 

<f>-=(w- k) + tp 2 Q £ l nT , 
with tp G Cq°(Qr), < (p < 1 such that (p = 1 on Q( Pl + P2 )/2, suppp C Q P2 and 

\D(p\ + (p 2 - pl)\pt\ < c/(p 2 ~ pi). 

Now, consider a cut-off function r\ S (Q( pi+p2 y 2 ) such that 0<r/<l, r/ = l 
on Q Pl and |D?7| < c/(p2 — pi). By Holder's inequality and the fact that \r]\ < 1 



\(w — /c)+| dxdt < c(n) / | (u> — /c) + 7y| dxdt 

Q P1 nn T ^S( P1 +p 2 )/2nn T 

< c(n) ( / \(w- A;)+?7| 2(1+ » ) dx dt J " +2 [F(p 2 , *)] 



2 

n + 2 



20 



P. BARONI, A. DI CASTRO, AND G. PALATUCCI 



To estimate the right-hand side we use Sobolev's embedding Lemma 2.11 and then 
twice Caccioppoli's inequality (3.19), 

f \(w- yfc)+r ? | 2(1+ " ) dxdt 

<c( sup / I (w — k)+(-, t) 1 2 da; J 

Vre(to-(pi+p 2 ) 2 /4,io) JB iPl+P2)/2 nn J 

2 1 f ,2 

\D(w — fc)+ dxdt+- / (u; — fc)+ dxdt 



: P2 nn T 1 V '" r| (P2 - Pi) 2 yg P2 nn T 



— c (t T9 / |(u; — /c) + | 2 dxdt 

V(P2 -Pir JQ P2 nn T 

we also used here the estimates for |r/| and |-Djy|. We estimate the first part of 
the latter term from above, keeping in mind that h < k and therefore A(k, p%) C 
A(h,p 2 ) 

I 1 2 / i 1 2 

|(tw — fe)+| dxdt < / |(w; — /i) + | dxdt 
P2 nn T JA{p 2 ,k) 

< / — h) + \ 2 dxdt. 

jQ P2 nn T 

Finally, it suffices to merge all the estimates above together with (3.18) (by replac- 
ing p\ with p2 there) to obtain the inequality in (3.17). 

At this point, we define the sequence of level sets and radii 

0<k i :=2(l-±)k u A:=| + 5 



2V 2 2 % 

where k\ is a positive constant which will be defined in few lines. We moreover 
choose into (3.17), for i £ IN fixed, p\ = Pi+i, P2 = Pi, ki = h < k = fcj+i. Call 
Ui := U(ki,pi) and Vi := V(ki,pi). By taking the a-th power, a > 0, of (3.17)2 
and the /3-th power, /3 > 0, of (3.17)i, and multiplying the resulting inequalities 
we finally have 

2 2(a+/3)i 

where c(a, /3) depends also on n, z^, L, Co, Po- Looking for i? > 1 such that a+/3 = ?9a, 
2a/ (n + 2) = i9/3 gives the value 



Q 1 1 2 



2 V 4 ra + 2' 
therefore, fixing a = (n + 2)t?/2 (and therefore /3 = 1), we have 



GLOBAL ESTIMATES FOR NONLINEAR PARABOLIC EQUATIONS 



21 



A well known iteration argument, see [19, Lemma 7.1], ensures that 
(3.20) lim JJfVi = if < ^HQj/fc?. 

i— >oo C 

Since, recalling that pi > R/2 and k{ < k±, 

0<([ \(w- k^+fdxdt) |Q fl/ 2 nO T n{w; > fci}| <U?Vi—>0; 

this yields w+ < k\ almost everywhere in Q.r/2 H Ot- We now find an appropriate 
choice of k\ ensuring that the condition in (3.20) is satisfied. Since 



I 2 1 2 11 

0-1= i/-r + 



4 n + 2 2 ra + 2# a 
recalling the definitions in (3.16) and the fact that p\ = R, we have 

/ r \a(#-l) 

[UiVl] = ( / Kw-^+fdxdt) \Q R nn T n{w > Ml*- 1 

\jQ.Rnn T J 

< \ QRnnrf~ 1 [ (w+) 2 dxdt 
jQ R nn T 

<\Qr\^ 1 f {w+fdxdt 



'Q R na T 

_J&a_ 

c 

if we choose as starting level 

■9 a 

2 C f 2 C f i2 

^ := — — — - / |w| dxdt = — — | / \w\ dxdt, 

\Qr\ JQ R nn T \Qr\ JQ R nn T 

c = c(n, L, Co, po). A similar argument, considering the negative part instead of 
the positive one, yields u>- < k\. Therefore we showed the validity of (3.15), but 
with the L 2 norm on the right-hand side. We show here how to lower this exponent, 
following [19, Chapter 7]. First notice that the estimate 



sup 

Q ri nn T 



i 

\w\ < c [ / \w\ dxdt 

\?2 — r\) n+l 



for all R < n < Ti < R holds as well; this can be seen with slight modifications 
of the previous proof, or adapting [19, Corollary 7.1] to our setting. Note that 
this last approach in our case would require a more involved argument: we should 
indeed distinguish the cases as in the first part of the proof of Theorem 3.5, and 
change a bit the radii. At this point, with this estimate at hand, we argue as 
follows: denote by U(r) the supremum of \w\ over Q r n Ot 5 for R/2 < r < R. Then 
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for R/2 < r\ < T2 < R and q > 0, using Young's inequality 

A well known iteration lemma (see, e. g., [19, Lemma 6.1]) at this point finally 
gives (3.15). □ 

We conclude this section by stating the counterpart of Lemma 3.2 up to the 
boundary. The proof plainly follows by that of Lemma 4.1 in [27] (see also [14, 
Lemma 3.3]), once taking into account the modifications needed in order to handle 
the fact that the cylinder intersects the boundary, see [40, Lemma 2.6] and the 
proof of previous Theorem 3.5. 

Lemma 3.7. Let u G C°([-T, 0]; L 2 (£l)) n L 2 (-T, 0; W 1,2 (O)) be the unique so- 
lution to the regularized problem (2.5), and let w be the unique solution of the 
Cauchy-Dirichlet problem (3.4). Then there exists a positive constant depending 
only on n,u such that 

/ (R-^u - w\ + \Du - Dw\) dxdt < cR / |/| dxdt. 



4. Proof of the main results 

This section is devoted to the proof of the main result, to the regularity analysis 
of the borderline case, and to precise Lorentz space estimates up to the boundary 
for the solution u. First, however, a couple of remarks allowing to simplify the 
proofs. 

Remark 4.1. Note that we can extend a solution to problem (1.1) to O2T := ^ x 
(— 2T, 0), where we also extend the function / to zero in £1 x (— 2T, — T] and the 
vector field a in the following way: a(x,t, z) := a(x, —2T — t, z) for t G (— 2T, —T). 
Therefore, denoting respectively by u, f the extension, we have that u satisfies 

Ut — diva(x, t, Du) = f(x, t) in f^T, 

U = on <9 par ^2T- 

Indeed, we briefly sketch this argument. Given a test function (p G Cq°(£12t), and 
supposing supp</?nfJ x {— T} 7^ in order to avoid trivialities, it is clearly enough 
to prove that the distributional formulation (2.2) is satisfied in C^(Qt) by the 
restriction of (p in $It, i.e. 

(4.1) — / u(ftdxdt+ / (a(x,t, Du), Dip) dxdt = / fipdxdt, 

Jn T Jq t Jn T 
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with ip as above. To this aim, take a C°°(]R) function ( e such that £ e = 1 on 
[— T + 2e, +00), ( £ = on (—00, — T + e] and \C' e \ < c/e, and put as a test function 
into (2.2) the product </> := (p( £ £ Cg°(Q T ). We get 

— / uiptCsdxdt— / u(/9^dxdi+ / (a(x,t, Du), Dip)Q e dx dt 
Jn T Jn T Jn T 

ftp( E dx dt. 



Now, we let e go to zero. Notice that both the first and the third integral on the 
left-hand side converge to the corresponding integrals in (4.1); the same happens 
for the integral on the right-hand side. It remains to consider the contribution 
involving the derivative of £ e ; for this, we have 



uipQ' £ dx dt 



-T+2e 



< c IMIl°°(o t )J _ / \u(x,t)\dxdt — 
thanks to (2.3). This completes the claim. 



This will allow to treat, by a simple argument, only cylinders not intersecting 
the initial boundary. Indeed for these, it would be enough to extend the solution 
as described above, and treat them as they were interior cylinders, or eventually 
cylinders just intersecting the lateral boundary. 

Remark 4.2. We stress here that we shall also not take into account the terminal 
boundary. Indeed the arguments in [5] could be also performed using not sym- 
metric, but backward parabolic cylinders as those in (2.1). This would apply in 
particular to parabolic cylinders of the form Q p (xo,0), with B p {xq) C £1. 

4.1. Integrability of Du. We firstly need a preliminary lemma, whose proof can 
be deduced by extending the original proof in the local case in [5, Lemma 6.2] (see, 
also, [40, Propositions 3.2 and 3.4]), which will permit us to apply the Calderon- 
Zygmund type Lemma 2.12 to the level sets of maximal operator functionals applied 
to the gradient of the solutions to our problem (1.1). Since the novelty here is given 
by the up to the boundary estimates, we prefer to state this result in the following 
form, that is in accordance with the analog in the elliptic case given by Proposition 
3.4 in [40]. 

Lemma 4.3. Let u be the solution to (2.5). Let A > 1 be an absolute constant 
depending only on n, L, v, cq, pq, and let xi > 1 be the higher integrability exponent 
as in Theorem 3.1 and 3.5 3 . Then for every S > 1 and any X > we have the 
following: if for some cylinder Q p , with p < mm{po, Ro}/12, po as in Definition 2.5 
and i?o the radius of a fixed ball Br , it holds 

\{(x,t) G Q p : M*(t nT \Du\)(x,t) > ASX}\ > S~ 2 ^\Q P \ , 



"A ctually, we take the minimum between the higher integrability exponent xi arising in both 
the aforementioned theorems. For the sake of simplicity, we keep the same symbol. 
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then there exists e = e(S,xi) such that 

Q P C {M*{1 Qt \Du\) > A} U {M*(l nT f) > eX} , 
where we denoted by M* := MX n and Mf := M? n . 

Proof. First, we notice that the statement can be restated as follows: 
Assume that for some cylinder Q p with p < min{po> ^o}/12 we have 

(4.2) Q p n {M*(t nT (\Du\) < A} n {M^l^f) < eX} + 0, 
then 

(4.3) \{(x,t) e Q p : M*(l nT \Du\)(x,t) > ASX}\ < S~ 2 ^\Q P \. 

By the hypothesis (4.2) there exists (yo, *o) 6 Qp = Qp(yi, si) such that, for any 
< r < .R , with (j/o, so) G Qr, 

(4.4) / ln T |Du|dyds < A 
and 

(4.5) r/ ln T |/|dyds<eA, 

reminding the definition of maximal operator. 

We want to prove that, for any (x, t) £ Q p (yi, si), the following inequality holds 

(4.6) M*(tn T \Du\)(x,t) < max{M*(l Q2p(n>sl)nQT \Du\)(x,t\ 3 N x}. 

For this, we consider a generic cylinder Q r such that (x,t) G Q r C Q_r ; for any 
f < P, Qr C Q2p(yi, s i)> and so we have 

-/ ln T p«|dyds < M*(l Q2p( , yijSl)nnT jZ)u|)(x,t). 

On the other hand, for any p < r < Rq we plainly have, enlarging, if necessary, the 
cylinder Q r , and using (4.4) 

/ ln T \Du\dyds < 3 N X. 

Thus, the inequality in (4.6) is proven. 

Now, let Qi p = Qi P (yi,si) D Q p and distinguish the following three cases: 
Q 4p C R n+1 \ fir, Qip n dp^T + and Q 4p C T . 

In the first case, the inequality in (4.3) does hold provided that A > 3 N , by 
taking into account (4.6). 

In order to deal with the second case, that is where the cylinder Q 4p intersects 
<9par^T> we extend to the parabolic setting the strategy in Proposition 3.2 in [40]. 
Precisely, take any (y2,s 2 ) G Sparer such that dist par ((yi, si), Sparer) := \yi - 
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J/2 1 + | si — S2 1 2 and consider the unique solution w to the following homogeneous 
Cauchy-Dirichlet problem 

{w t - div a(x, t, Dw) =0 in Qi 2p (y 2 , s 2 ) n 
w = u on 3 par (Qi2 P (y2, s 2 ) n n T ). 

Also, we extend by zero the function u to \ Q,? and by u the function w to 

H" +1 \(Qi2p(y2,s 2 )nO r ). 

Now, we are in a position to estimate the measure of the investigated level sets 
by using (2.9) twice with 7 = 2%i and 7 = 1. It follows 

| {(*,*) € Q p : M*(l Q2p{yusl)nUT \Du\)(x,t) > ASX}\ 

< \{(x,t) € Q p : M*(l S2p(j/liSl)nQT |^|)(x,t) > ASA/2} I 

+ |{(x,t) € Q p : M*(l Q2p[yi>Sl)nnT \Du- Dw\)(x,t) > ASX/2}\ 

(4.7) <c(ASA)" 2xi / \Dw\ 2xi dyds 

jQ2 P (yi,si)nn T 

+ c(ASA)~ 1 f \Du- Dw\dyds. 

By Theorem 3.5, applied with R = Yip < po, the first term in the right-hand side 
of (4.7) is bounded from above by 

c(AS\)- 2 ^\Q 6p \(—}—[ \Dw\dyds) , 

\jyi2p| J Q 12p (y 2 ,s 2 )nn T J 

since Q 2p (yi,s 1 ) C Q6p(?/2,s 2 ). Thus we get 

\{(x,t) 6 Q p : M*(l Q2p{yuSl)mT \Du\)(x,t) > AS\}\ 



< c(ASX)- 2 ^ \Q 6p \l I \Du\ dy ds ) 

\\<dl2p\ JQ 12 p(y 2 ,s 2 )nQ T 

+ c{ASX)- 2xi \Q 6p \l-^—[ \Du-Dw\dyds ) 

\ \Ql2p\ JQ 12p (y 2 ,s 2 )nn T J 

+ c(ASA)~ 1 ^4 / \Du- Dw\dyds. 

\Ql2 P \ JQ 12 p(y 2 ,s 2 )nn T 

In view of the comparison estimates established in Lemma 3.7 (taking R = Yip < 
po there), together with (4.4) and (4.5), since (y ,s ) £ Q P (yi,si) C Qi2 P {y2,s 2 ) C 
Qrq, and choosing e = S _2xi+1 , we arrive at 

\{(x,t) e Q p : M*(l Q2p{yuSl)nnT \Du\) > ASA} I < c(A" 2 ^ + A" 1 )^^ \Q P \. 



\ 2X1 
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Now, it suffices to choose the constant A so that A > 3^ and 2cA 1 < 1/2; i. e., 
A > max{3 N , 4c}. It follows 

|{0M) € Q P : M*(t Q2p{yuSl)nnT \Du\)(x,t) > ASX}\ < V^jQpl, 

which, again in view of (4.6), yields (4.3). 

Finally, when Q^ p C Qt, the inequality in (4.3) plainly follows by Lemma 6.2 
in [5], hence using the local interior results given by Theorem 3.1 and Lemma 3.2. 

□ 

Before starting with the proof of Theorem 1.2, we want to emphasize that most 
of the differences between the problem we are dealing with and the analog in the 
elliptic case analyzed in [40] have arisen in the previous section. For this, the 
general strategy of the proof follows that of the proof of Theorem 6.1 in [5] (see, 
also, [14, Theorem 4.1]. We prefer to give some details for the reader's convenience, 
and we sketch the proof in three steps. 

Proof of Theorem 1.2. We denote by Qo = Bh (xq) x (to — Rq, to) a fixed cylinder 
in R n+1 that contains Qt, with its basis Bq = Brq being a ball of radius Rq < 
2diam(f2). Note that since Qo D we can assume Rq > diam(fi)/2. Now we 
consider the following maximal operators 

M*(- ) := M * Qo (lfi T • ) and M*( ■ ) := M* t Q Q (lfi T • ). 

Step 1 - Application of C alder on- Zygmund-Krylov-Safonov covering type lemma. 
We consider the following sets X and 3^ defined 

X := {(x,t)eQ :M*(\Du\)(x,t)>(AS) k+1 \ 1 }, 

y := [(x,t) € Qo : M*(\Du\)(x,t) > (A5) fc Ai 

and M*(f)(x,t) > e(AS) k X i y 

with 

(4.8) Ai := cS 2xi j^-r [ \Du\dxdt, 

\Qr\ JQ nn T 

where S > 1, A, xi, £ are as in Lemma 4.3. We want to apply Lemma 2.12 with R = 
minjpo, Rq}/12. Firstly, we notice that by Lemma 4.3 (choosing A = (^45) fc Ai > 
there) the assumption (ii) is satisfied with 5 = S~ 2xi . Thus, it remains to prove 
that \X\ < 5 _2xi |Q^|. For this, it suffices to use the boundedness of the maximal 
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operators in Marcinkiewicz spaces together with the choice of Ai in (4.8). Indeed, 
{(M) 6 Qo : M*(\Du\)(x,t) > (AS)* +1 Ai} 

\Du\ dx dt 



< 



<S^\Q R \, 

where we also used the fact that A, S > 1. Hence, the application of Lemma 2.12 
and the definitions of X and y yield 

{OM) e Q : M*(\Du\)(x,t) > (AS) fe+1 Ai} 

< c 3 5- 2x \{(x,t) G Qo : M*(\Du\)(x,t) > (AS) k X x } 

+C3 \{(x,t) G Qo : M*(f)(x, i) > e(AS) fc Ai} , 
for every G N; with c% depending only on n, po, and diam(fi). 

(fc + l)jV 7 

Now, by taking into account a multiplication by a factor (AS) N ~"i , from the 
previous estimate we deduce 



(fc + l)iV-y N ~t , , 



(4.9) 



N-yk N'y N-y 



< c 3 (AS)^ A^ S^~ 2x \?- J ^((ASfXt) 

N-yk ( AS\ N ~t 1 1 

+ c 3 (AS)— - (X^—^AS)^!), 



where, for any K > 0, we denoted by 

(4.10) m(K) := \{(x,t) G Qo : M*(|£>u|)(sM) > #}| 
and 

(4.11) /. 2 (iO := | G Qo : M?(f)(x,t) > K}\. 

Now, we note that, since xi > 1 and 7 < 2N/(N + 2), we have that the quantity 

iV 7 



d := 2 X i 

TV — 7 

is positive. Therefore, we can choose 



> 2(xi - 1) > 



JV7 

S := [AA^ 



and, by computations, inequality (4.9) provides the existence of a constant c = 
c(n, L, v, diam(fi), po) such that, for every > 0, 

(4.12) (AS 1 ) - ^ A 1 w_ ' , /ii((i45) fc+1 A 1 ) 



N-yk 



N-yk 



jV 7 



< ^(AS)^X^ ^((ASyXx) + c(AS)^(eAi)^^ 2 ((^5) /c eAi). 
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Step 2 - Level sets estimates. In order to establish some Lorentz spaces estimates 
on level sets we proceed as follows. We take < (3 < oo and operate some compu- 
tations to get 

' J^L. /3(iV-7) HA 

[X N -y /j,i(X)] N i — 
A 

< Q + cP \ASf log(AS)) X{\QdT^ + ~c?{ASf \og(AS) J {oo), 



where 



J(oo) := {(AS^ieXx^^AS^eXx) 

k=0 

and the constant c > 1 is increasing in the variables n, L, v and decreasing in f3, 
such that c — > oo as /3 — > 0, while it remains bounded when /3 is bounded away 
from zero. Thus, from this inequality, we can plainly deduce 

J [X—^(X)]-^T^ < l- + 2c (3 (AS^log(AS)\X^\Q \^T- 

(4.13) +cP{ASfP / [A— ^(A)]^-^. 

Jo A 

We recall the definitions of \i\ and fi2 given in (4.10)-(4.11), and, now, we choose 4 
/3 = q G (0, oo), so that, by the definition of parabolic maximal operators, we get 

AT 1 

(4.14) \\M*(\Du\)\\ Li ^ q)iQo) < c^lQol - + c||M 1 *(/)|| Lte)g)(So) , 

up to relabeling the constant c, by keeping the same properties as before. 

Now, thanks to a classical global estimate for the gradient of u established in [9] , 
which we will use as in the form presented in [4, Lemma 5.1] where the exact 
dependence on the L 1 -norm appears, we can estimate Ai from above, as follows 

Ai < c(n,p ,diam(0))|Qo|^ _1 / \f\dxdt. 



N-/ 



T 



Therefore, for any (x,t) € Qq we have Mf(f)(x,t) > Xi/c(n, po, diam(f2)), and 
thus 

jv — i 

A,ia,i - < ll^(/)ll l(A)(fin) - 

Finally, in view of the standard inequalities \Du(x,t)\ < M*(\Du\)(x, t), for a. e. 
(x,t) e So, and |[Mf (/)|[ £ / ^ \ (£o) < ||/||i( 7 , g )(Q ), we arrive at 

( 4 - 15 ) \\ Du \\ L ^ >q yn T ) ^ c ll/ll^(7, 9 )(So)> 



4 We note that in the computations above we preferred to keep the generic exponent /3, since 
it could be chosen arbitrarily in order to obtain estimates also in other rearrangement and non 
rearrangement spaces. 
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where the constant c depends only on n, L, v, 7, q, po and diam(Q). Recalling that 
r^T C Qo an d / = in R, n+1 \ Qt we plainly deduce the desired estimate in (1.4). 

Similarly, we can deal with the Marcinkiewicz case q = 00 (see, e.g., Step 3 in 
the proof of [5, Theorem 3]) and we arrive at 



\\Du\\ Nj_ < C\\J ||A47(f2 T ). 

Step 3 - Conclusion of the proof. We recall that we proved the estimate in (4.15) for 
the approximating solutions u = u/- to problem (2.5) with / = fk = ?&(/), where 
Tfc is the truncation operator defined by (2.6). In order to conclude, it suffices 
to use the lower semicontinuity of the Lorentz norms together with the standard 
approximating arguments stated in Section 2.1. Hence, we have proven that 

( 4 - 16 ) \\ Du k\\ L f_N2_ jq \ (nT) <c||/k|U( 7 ,<r)(Co), 

where the constant c does not depend on k. We also notice that in (4.16) we have 
used the fact that 



Wfk\\L( 7 ,q)(Q ) < H/llL(7,g)(Qo)' 
since by the definition of fk it holds that \fk\ < [/[. Thus, in order to pass to the 
limit on k — > 00 in (4.16), it suffices to use the lower semicontinuity of the Lorentz 
norms (recall remark on page 8) together with the approximating arguments stated 
in Section 2.1 (recall, in particular, (2.6) there). □ 

Remark 4.4. As in the classic case, Theorem 1.2 fails for the borderline choice 7=1. 
Even in the elliptic case, one has to impose some further LlogL integrability on 
the datum / in order to obtain the following implication 

1 N 

(4.17) / £ L (— T, 0;LlogL($7)) \Du\ £ L~ (ft T ); 

see Theorem 1.8 in [8]. However, we could extend the potential approach in the 
proof of Theorem 1.2 to the analysis of the borderline case in the parabolic Orlicz 
space LlogL(Oy), to get, for any cylinder Q p C R n+1 , the following estimate 



J_ JV-l 
AT II #11 JV 



(4.18) '^""ii^TcCp/anftr) - "Li(Q P nn T y\J i\Li og L(Q p nn T y 

We stress that when considering the preceding estimate in the whole fif> the by- 
now classic result (4.17) in [8] is stronger, since it does not require any regularity 
assumptions on the domain -T2 and it does hold in the weaker assumption of / being 
in L l (—T, 0; L log L(U)), despite (4.18) provides local estimates up to the boundary. 



4.2. Integrability of u. This section is devoted to Lorentz space estimates for 
the solution it to (1.1). We will use the same techniques used to prove the spatial 
regularity of the gradient. Hence, we will obtain an estimate on the level sets of the 
maximal operator associated to u, in terms of the level sets of a maximal operator 
of the datum /, up to a correction term; that is, the equivalent to (4.14) for u. This 
will allow to deduce the proof of Theorem 1.4, stated in the introduction, which is 



30 



P. BARONI, A. DI CASTRO, AND G. PALATUCCI 



close to the one of Theorem 1.2, but also a higher regularity result is needed, as 
stated in the second part of Theorem 3.5 proven at the end of Section 3. 

Keeping in mind the notation used in Theorem 1.2, we have the analog of 
Lemma 4.3 for u. Again, we recall that in what is following we always deal with 
the approximating solutions u = Uk defined in Section 2.1, and with the datum 

/ = fk- 

Lemma 4.5. Let u € C°([-T, 0]; L 2 (0)) n L 2 (-T, 0; W 1,2 (n)) be the solution 
to (2.5). Let A > 1 be an absolute constant depending only on n,L,is, and let 
X > 1 be the higher integrability exponent as in Theorems 3.1 and 3.5. Then for 
every S > 1 and any A > we have that exists e = e(S, x) such that the following 
statement holds: If for some cylinder Q p with p < min{po, Rq}/12 it holds 

\{(x,t) € Q p : M*(u) > ASX}\ > S~ 2 *\Q P \ 

then 

Q p c{M*(u)>\}U{MZ(f)>£\}, 
where we denoted by M*(-) := M£ Qr (ln T -) and M|(-) := M 2 * Qr (ln T -)- 

Proof. The proof of this lemma is essentially the same as that of Lemma 4.3. It 
first suffices to substitute M*(\Du\) and M*(f) by M*(u) and M|(/), respectively. 
Then, we need to use the comparison estimates for u stated in Lemma 3.2 and 
Lemma 3.7, and the higher integrability, again for u, given in (3.3) and (3.6). This 
will permit to introduce the parameter x > 1 which can be chosen arbitrarily large; 
and this is another difference with respect to the estimates in Lemma 4.3, in which 
such a parameter (read xi > 1) was fixed. □ 



Proof of Theorem 1.4- Since we closely follow the proof of Theorem 1.2, we prefer 
to keep the division of the proof in the same three steps. 

Step 1 - Application of Calderdn-Zygmund-Krylov-Safonov covering type lemma. 
We want to apply Lemma 2.12 to the sets 

X := {(x,t)eQ :M*(u)(x,t)>(AS) k+1 X }, 

y := {(x,t)eQ :M*(u)(x,t)>(AS) k \ 

and M|(/)(M) >e(AS) k X }, 

with A :=cS 2 *-L f \ u \ dxdt and ^ = min{p ,fi } 
\Qr\ Joanna 12 



In view of Lemma 4.5, the assumption (ii) of Lemma 2.12 (choosing A = (^4S') fc Ao 
there) is satisfied with p < R = min{poi ^o}/12 and 5 = S~ 2x . To prove (i), that 
is \X\ < aS' _2x |Qo|) it suffices to use the boundedness of the maximal operators in 
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Marcinkiewicz spaces together with the choice of Ao and the fact that A, S > 1. It 
follows 

iV 7 (fc + l) N ~/ , , 

(AS)^^X^^((AS) k+1 X ) 

Njk Ni N-y _Ny_ 

(4.19) < (AST)N-*AN-*S»-*~ 2xo \g- ar 'iii((AS) k \Q) 

N-y N-yk N-y . 

+ (AS/e)^^ {AS)~ (e\ )~ n 2 ((AS) k eX ), for every k > 0, 
where, for any K > 0, we denoted by 

IH(K) := \{(x,t) € Q : M*(n)(x,t) > K}\ 

and 

// 2 (i<0 := |{(x,i) e Q : M 2 *(/)(z,i) > K}\. 

At this level, we can take advantage of the possibility to choose x large enough to 
satisfy 

iV 7 

d := 2 X - — — — > 0, 

/ j N-y N-y r, 

and, using the 5 := UA^- 2 ^ , it follows A N ~ 2 i S N ~ 2 ~< x < 1/4, so that the 



estimate in (4.19) becomes 

Nj(k+1) N ^ , , 1 N-yk N ~> 

(AS)^r X£- 2 ~< vi((AS) k+1 X ) < -(AS^X^^dASfXo) 

N-yk N-y , 

+ c(AS) — (eX ) — fi 2 ({AS) k eX ), 

for every k > 0. 

.Step 2 - Level sets estimates. Arguing as in Step 2 in the proof of Theorem 1.2 we 
obtain the following Lorentz estimates for the solution u, 



for any < q < 00. Similarly, in the case of Marcinkiewicz space, that is q = 00, 
we arrive at 



Notice that at this stage we used the fact that M| has a higher regularizing effect 
with respect to that of M£. Moreover since we are assuming / = in R n+1 \ Q,t 
and Qt C Qo we obtain exactly the desired estimate. 



Step 3 - Conclusion of the proof. We can now conclude the proof using standard 
approximation methods, again by means of the lower semicontinuity of the Lorentz 
norms. □ 
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